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• which leads to a correct exponential decrease of the amplitude at large impact parameters, affects the main 

(~\ properties of the BFKL Pomeron. We solve the equations for total cross section and (|6 2 |) numerically and 

C*~> developed some methods of analytical solutions. The main result is that the modified BFKL Pomeron has 

. the same intercept and a'jp = as the BFKL Pomeron. We found the energy dependence of the saturation 

. _ scale which turns out to be much milder than for the BFKL equation. 

X: 



Keywords: BFKL Pomeron, solutions to the BFKL equation, semi-classical approach, diffusion 



approximation, large impact parameter behaviour of the amplitude . 



PACS: 12.38-t, 12.38. Cy,l 2.38. Lg, 13.60.Hd, 24.85. +p, 25.30.Hm 



Contents 

[I]. Introduction [I] 

2. Impact parameter dependence of the BFKL Pomeron |^ 

p7l| The BFKL Pomeron: generalities y 

U Equation for (|6 2 (Y, /) |) | 

p^ Numerical solution || 

|. Modified BFKL Pomeron | 

|Q| Pomeron intercept [To| 

[3.1. 1| Numerical solution for the Pomeron inrtercept [l0| 

[3.1.2| Variational method [0] 



3.1.3 Semi-classical solution [U3] 

3.1.4 Diffusion approximation [l5| 
Pomeron slope |l7] 
Saturation momentum O 



Conclusions 20 



|||. Acknowledgements 21 



1. Introduction 

The large impact parameter dependence of the scattering amplitude has been the principle but still un- 
solved problem in the CGC/saturation approach for the past decade. It was shown in Refs. [1-4] that 
CGC/saturation approach [5-8] that leads to the partial amplitude smaller than unity and satisfies the 
unitarity constraints, generates the radius of interaction that increases as a power of energy in explicit 
contradiction to the Froissart theorem [9]. It stems from large b behaviour of the BFKL Pomeron [10,11] 
which has the form: A(b^> 1/Q S ) oc s /b 2 *. Amplitude A(b^> 1/Q S ) becomes of the order of unity 



*The more detailed discussion of the impact parameter behaviour of the BFKL Pomeron will be done in the next section. 



at typical b 2 oc s leading to a oc s in the contradiction to the Froissart theorem (<r oc In s). The 
power-like dependence of the scattering amplitude is a direct consequence of the perturbative QCD tech- 
nique which is a part of the CGC/saturation approach. Since the lightest hadron (pion) has a finite mass 
(to?,-) we know that the amplitude is proportional to exp (—2m n b) at large b. This exponential behaviour 
translates into Froissart theorem. Therefore, we have to find how confinement of quarks and gluons being 
of non-perturbative nature, will change the large b behaviour of the scattering amplitude in the region 
where this amplitude is small. 

This complicated problem in spite of numerous attempts [4, 14-20], has not been solved. However 
we learned several lessons from these tries. First, in the framework of the DGLAP equation [21] we can 
factorize out the non-perturbative large b behaviour writing for the scattering amplitude^ 

A (b, Y,r) = S (b) A DGLAP (Y, r) (1.1) 

(see Ref. [22] ). Indeed, considering the scattering amplitude at fixed transferred momentum q ( which is 
Fourier conjugated to b), one can see that for q < /J. S oft the evolutions in ln(l/r) do not depend on q. 
However, for q > /J, so ft the logs take the form ln(l/(rg)) and the q dependence cannot be absorbed in 
S (b) in Eq. ( |1.1[ ) [22]. Using Eq. ( |1 . 1[ ) we can absorbed the non-perturbative corrections at large b in the 
definition of the saturation scale Q s (Y; b) [18, 23-27]. 

However, such way of including the non-perturbative large b behaviour does not work [4, 14-17] in the 
case of the BFKL and BK evolutions [11,28]. Since we are interested in the behaviour of the scattering 
amplitude at large b where this amplitude is small, we need to find a way to introduce the non-perturbative 
corrections directly to the BFKL kernel. Hence the non-linear dynamics does not influence on a solution 
to this particular problem. We would like to recall that the saturation scale and its dependence on b 
follows directly from the solution of the BFKL equation (see Ref. [8] and reference therein). It has been 
check by numerical calculations (see Refs. [14-17]) that if we modify the BFKL kernel introducing by 
hand a function that suppressed the production of the dipoles with sizes larger than l//j, so f t , the resulting 
scattering amplitude has the exponential decrease at large impact parameters. 

In this paper we modify the BFKL kernel in the following way: 
a s K BFKL (x 13 ,x 32 \x 12 ) = «5^V => a s -^e- B (*h + *h) = a s K B (x 13 ,x 32 \x 12 ) (1.2) 

x 13 x 32 x 13 x 32 

Having this modified kernel we are going to answer the following questions :(i) how the intercept of 
the BFKL Pomeron depends on B; (ii) what is the dependence of (|6 2 |) on Y and the size of the dipole; 
and (iii) what is dependance of the residue of the BFKL Pomeron on the size of dipole. The goal of this 
paper to compare the BFKL Pomeron with the modified kernel to the soft Pomeron we know both from the 
Regge high energy phenomenology [29,30] and N=4 SYM with AdS-CFT correspondence [31-33]. These 



'In this paper we use the following notation: Y = ln(l/:r) where x is the fraction of the energy carried by the dipole, 
r is the size of scattering dipole, q is the momentum transferred for the scattering amplitude and fJ. ao ft is the scale of soft 
interaction (fj, so ft ~ ^-qcd)- Notice that q is the Fourier conjugated to the impact parameter b. 



approaches leads to the soft Pomeron with sufficient large values of the intercept and with the slope (a'jp) 
which is equal to zero ( (|6 2 |) oc a'jpY). 

The result of the paper are the answers to these three questions. We found out that the intercept 
for the modified BFKL Pomeron is the same as the intercept of the BFKL Pomeron with original kernel 
(B = in Eq. fll.2j ))i . At high energies (\b 2 \) — > Const. In other words we expect that a'jp — > at large 
Y. The Pomeron residue does not depend on the dipole sizes (r) for r < 1/B but it drops for r > 1/B. 
In short we see that the Pomeron with the modified kernel matches the soft Pomeron as we know it both 
from N=4 SYM and phenomenology. 

The paper is organized as follows. In the next section we consider the BFKL Pomeron and discuss 
its main properties concentrating our attention mostly on the impact parameter dependence. In section 
3 we present the numerical solution for the modified BFKL Pomeron with the kernel of Eq. (|1.2| ) with 
B ^ 0. In this section we develop several analytical methods to evaluate the intercept of the modified 
BFKL Pomeron: variational method, semi-classical and diffusion approximations. We solve the equation 
for (\b 2 \) and show that the numerical solution and the analytical estimates lead to (|6 2 |) which does not 
depend on energy. In addition, we evaluate the saturation momentum which turns out to show much milder 
energy behaviour for the modified BFKL Pomeron than for the BFKL equation. 

In conclusions we summarize the results and compare with the soft Pomeron. 

2. Impact parameter dependence of the BFKL Pomeron 

2.1 The BFKL Pomeron: generalities 

The general solution to the BFKL equation for the scattering amplitude of two dipoles with the sizes n 
and T2 has been derived in Ref. [10] and it takes the form 

N(r u r 2 ;Y,b) = (2.1) 

E f^ L -4n(T,r2) d 2 R Y d 2 R 2 d(R 1 -R2-b)e^^ Y E^' n (r 1 ,R 1 ) i? 1 ^'" (r 2 ,i? 2 ) 



n=0 



with 



w( 7 ,n) = a sX (7,n) = a s (2^ (1) - V> (7 + M/2) - V (1 - 7 + M/2)) ; (2.2) 

where ^ (7) = dlnT (j) /dj and T (7) is Euler gamma function. Functions E n ^ (pi a ,p2a) are given by 
the following equations. 

(\ 1— 7+n/2 / * \ 1— 7— n/2 
-^-) -^r) , (2.3) 

PlaP2a) \p\aPla) ^ ' 

In Eq. ( |2.3| ) we use the complex numbers to characterize the point on the plane 

Pi = %i,i +ixiX, P* = Xi,i~ixi,2 (2.4) 



'We will call the BFKL Pomeron the solution to the equation with the kernel of Eq. (IL4) with B = 0. 



where the indices 1 and 2 denote two transverse axes. Notice that 

P12P*12 = rf; PlapXa = [&i ~ ifi) P2a P* 2a = 1^4 + ^i) (2-5) 

At large values of Y the main contribution stems from the first term with n = 0. For this term Eq. (|2 
can be re-written in the form 

E^ (ri,Ri) = (— ^-_ j . (2.6) 



The integrals over R\ and R 2 were taken in Refs. [10, 12] and at n = we have 

H^(w,w*) = |d 2 J2i^' (ri ) fli)>£ 1 ~ 7,0 (r2,£i - &) = (2-7) 

( 7 (1 _ 2 7 ))2 { 6 7 ^ 7 ™*~< F (7, 7, 2 7 , w) F (7,7,27,™*) + 
&i_ 7 u; 1 -W 1 ~ 7 .F (1 - 7, 1 - 7, 2 - 27, to) F (1 - 7, 1 - 7, 2 - 2 7 , u>*) j 
where F is hypergeometric function [13]. In Eq. ( |2.7| ) wtf* is equal to 



WW = o 9 ^" 



^-i(fi - f 2 )) (6+ i(ri - r 2 )) 
and 6 7 is equal to 



r( 7 ) r(i- 7 ) 
( 

Finally, the solution at large Y takes the form 

(2.10) 

Eq. ( |2.10| ) shows that at large b S> T\ and r 2 ww* = r\r\jb <C 1. Therefore, we can replace F 



b _ ^3 2 4(l/2-7) U7J M1-7J f2 gs 

7 r(i/2-7)r(i/2 + 7 )- (2 - 9j 



Ar(ri,r 2 ;y,6) = / ^ 4n (r, r 2 ) e w(7 ' 0) Y # 7 (w, w*) 

I 2 IT l 



functions in Eq. ( [2.7] ) by unity and Eq. ( 2.10 ) degenerates to the following expression 
nn,r 2 ;Y, b) = /^«(7;, 2 )^-{6 7 (^) 7 + 6 1 _ 7 (^) 1 " 7 } ^ 



ni2_^ y (211) 



b 2 

where at F > 1 7 ->■ 1/2 and w = a5x(l/ 2 )- One can see that AT (ri, r 2 ; Y, b) < 1 for 6 2 < rir 2 e woY [1-3]. 
For DGLAP evolution the essential 7*1 <C r 2 and 7 — )• and Eq. ( 2. 10] ) takes the form 



N( ri ,r 2 ;Y,b) = [ ^L <^( 7 ;r 2 ) e^ Y (ww*)" (2-12) 

J 2TTI 

One can see that or b <C r 2 w w* = r\jr\ and the impact parameter dependence can be introduced 
through non-perturbative initial condition. However, for r\ <C \b — 2^*2) "C r 2 ww* = r\/\b— t^I 2 and 
b dependence cannot be absorbed in the initial condition. 



2.2 Equation for {\b 2 (Y,l) |) 

In this section we derive the equation for (|6 2 (Y, I) |) defined as 



(\b 2 (Yl)\) - fdWN—Hn,rr,Y,b) 



The BFKL equation takes the form: 
dNBFKL ^ b] Y) =a s j d 2 x 13 -*kr { 2 NBFKL (**. & - ^32; y) " iV B ^ L (x 12 ,6;y)} 

J x 32 1 \ l J x 13 > 

where N BFKL (212, 6; Y) = N BFKL (212, 6; Y) /x 2 2 and X12 is the size of the dipole (n in the notation of 
the previous section). The size of the second scattered dipole r 2 we suppress in the notation. 

Integrating Eq. ( p.!4| ) over the impact parameter we obtain the equation for J\f BFKL (x 12 ; Y) = 
/ d?bN BFKL (x 12 , b; Y) which takes the form 



dM BFKL (x 12 ;Y) 



as [ d 2 x 13 4" {2AA B ^ (z 13 ; Y) - ^W*™ (*125 ^ } (2-15) 

J X 32 L x 13 J 



oY ,, 0.32 <- j, 13 



Multiplying Eq. (J2.14Q by 6 2 we derive the equation for J\f BFKL ( Xl2 ; Y) = J d 2 bb 2 N BFKL (x 12 , b; Y): 



<9Y 



(2.16) 



= a s // dH> d 2 x 13 -i-^+k 2 ) 2 4- (2 iV fl ™fxi3, V^b-lx^Y) - X -^N BFKL {x 12 ,b;Y)\ 

= a 5 /d 2 x 13 ^{2^^(x 13 ;Y) - ^™^ 12 ;Y)} + ^a s [ d 2 x 13 Af BFKL (z 13 ; Y) 
J x 32 k x 13 J I J 

+ U«s f[d 2 b'd 2 x 13 9-x 32 -i-2 N BFKL (x l3 , b' = b-\x 32 ;Y) =ol (2.17) 



X 32 



The last term (see Eq. ( 2.17 ) ) is equal to 0. It follows directly from Eq. (|2.10| ) and the expression for 
ww* of Eq. fl2JD, since they show that N BFKL (xi 3 ,9; y) is even function of 9 (N BFKL (x 13 ,9;Y^ 



jyBFKL l Xl3j —b';Y)). Therefore, the integral over b' in Eq. ( |2.17j ) vanishes. 



We need to solve Eq. ( gig) and Eq. (|I|) to find {\b 2 (Y, l)\) = M BFKL (x 12 ; Y) M BFKL (x 12 ; Y). 



The initial conditions for these equations are taken in the form 

M BFKL (x 12 ; Y = 0) = In (1/ (x 2 12 A 2 QCD )) ; M BFKL (x 12 ; Y = 0) = b 2 In (1/ (x? 2 A| CD )) ; (2.18) 

where 6q is the value of (|6 2 |) at Y = 0. 



Eq. ( 2,15 ) has a well known solution: 



e+joo 



dv 



N BFKL (x l2 -Y) = / 5p L_ e w M y + (-M' (2.19) 

where/ = In f xf 2 Aq CZ) J and w (y) = w(|+w, 0) (see Eq. (|2~2D). 

At large y the main contribution stems from i/ — ^ where we can use the simplified form of BFKL 
kernel u (7,0): its expansion at small values of v (diffusion approximation) 

u{y) = wo - D Q v 2 with uj = 41n2a 5 = 2.772 a S and D = U( (3) a s = 16.828 a s (2.20) 
Taking integral over v we have 

*iff KL (xn;Y) = -^=e^ Y -^-^ (2.21) 



The second term in Eq. (|2.16j ) is a function of only Y . We solve Eq. ( |2.16| ) using double Mellin 
transform 

f-e+ioo i r-ie+ca 

N BFKL {x 12 -Y) = / pL n(co,v)e" Y + (~^) 1 (2.22) 

Plugging Eq. fl2.22|) in Eq. (p. 16 ) we obtain that 



( u- uM )*M = Wi/ 1' „ ', /<«^+~'>< 



2 



«/) 7 2vrio; - w (i/) (!-«/')' 



I _ it b% f dv' I 1 1 

2 aS (p^) 7 2^-^') (iz/-±) 2 JT^X (2 ' 23) 



The l.h.s. of Eq. ( 2.2$ ) is the result of the integration over X13 in Eq. fl2,16|) . This integration demonstrates 



that in the framework of the BFKL equation we cannot find j\[ BFKL (x\2~, Y) since the integral is divergent. 
It is expected since the BFKL equation is conformal symmetric and therefore, the dimensional scale can 
be originated only in the initial condition. For example, in dipole-dipole scattering the natural scale is the 
size of the target dipole (bo)- In this case, the choice of I is I = In (xf 2 /^o) w hich is used in Eq. ( [2.23|) and 
the initial condition of Eq. Q2.18 ) can be written as 



M BFKL (x 12 ; y = 0) = b 2 In (6g/ (z? 2 ) ) ; (2.24) 

As one can see in Eq. ( |2.23| ) we integrated over I < and, thereby, introduced the infrared cutoff, 
assuming that the dipole size is less than bo . Hence Eq. ( |2.23| ) includes a modification of the BFKL 
equation, which we apply only for the dipoles with sizes smaller than 60 • For the modified BFKL equation 
such a cutoff is an intrinsic property of the equation. Therefore, introducing the cutoff in Eq. ( 2.23J ) we 



model the modified BFKL equation with the kernel of Eq. (1.2) with the goal to obtain the analytical 
solution which we cannot hope to get for the modified kernel of Eq. ( |1.2| ). It should be noticed that the 
factor 1/(2 — iv) stems from the condition that the non homogeneous term in Eq. ( [2.23 ) we consider only 
for I < 0. 



Generally speaking we can add to solution of Eq. ( 2.23 ) any solution of the homogeneous equation in 
the form 



An (u>, v) = n (v) 



1 



(2.25) 



UJ — UJ (u) 

with arbitrary n (y). 

uj (v) has pole of the first order at iv —> ±1/2 and at large values of Y the main contribution stems 
from v — > 0. Using this features of uj (v) and its expansion at small values of v (see Eq. ( |2.20| )) we can take 
the integral over v' in Eq. ( |2.23| ). It takes the form 



1 

-7ra 5 — 



bl 



a s 



-,-iv) \as \/Dq (uj - uj ) 
Using Eq. ( |2.26 ) the general solution can be written as 

I b 2 I / vr 2-Ka.s 



n (uj, v) 



+ 



uj - u(y) (\-iv) \ \2 ' y/Do (uj - uj ) 
Function n (v) should be found from the initial condition 



n (u) 



e+ioo 



duj 



7T 



oo 2ni uj — uj (u) 
Resolving Eq. ( p. 28 ) we obtain 



+ 



2iras 



7a^ 



UJ - UJq) 



n (v) 



iv 



n (u) 



7T 

2 



iv 



2iras 
~D^v~ 



Finally 



n (uj, v) 



ti 2 



uj - uj(v) (i -iv) 
Taking integral over uj we get for v — > 



+ 



2iras 



€+io ° duj _ , 
ir-. n{ ,u,v) 



h 2 



M»)Y 



2 ~ iv V D (W - WO) 

1 2-7TQS \ 2(Xg', 



2na s \ 
D v J 



iv) 



+ 



?vr 



iv 



D v J 



DnV 



Erfc VD^Y 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



One can see that at small Dqv 2 Y the factor in curly bracket vanishes. Taking the integral over v using 
the steepest decent method we see that this parameter is equal to D^v^pY = l 2 / (4JDqY). Since we are 
interested mostly in the solution at I which is not very large concentrating our effort on so called Regge 



domain, we can safely consider this parameter as being small. Expanding Erfc {^/Dqv 2 Y^ with respect to 
this parameters, we see that < . . . > = 1 + (2asy/Tr / \/Do) VY. 

Therefore, plugging solution of Eq. ( p. 31 ) into Eq. (|2.13| ) one can see that 



(\b 2 \) 



% 1 + 



fJT n 



(2.32) 



in this approximation. Such a behaviour of (\b 2 \) versus Y has been expected (see Ref. [39]). Comparing 
this behaviour with (|6 2 |) = Aa'jpY one can see that a'p for the considered modification of the BFKL 
Pomeron is equal to zero. 
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Figure 1: d In Uf BFKL (I- Y) \ jdl (see Fig. @-a) and In Uf BFKL {l\Y)\ (see Fig. |-b) for the BFKL equations (see 
Eq. ( |2.15D and Eq. ( |2.16 ) with the BFKL kernels) as function of I at Y = 43. as is chosen to be equal to 0.2 and 
/ = In (s^ 2 Aq CD ). Curves correspond to different values of R. The red curve shows the solution of Eq. ( 2.21 ). 



2.3 Numerical solution 



Searching for numerical solution of Eq. ( 2.15| ) we introduce the regulator at short distances R in the 
following way 



dM 



BFKL 



{xi2-,Y) 



dY 



as d xi3 



.32 



+ .R 2 



2 AT 



BFKL 



x 12 



ln:CV, - > , , ::\ n R2 N BFKL (x 12 ;Y)} 

(2.33) 
We solve Eq. fl2.33|) at fixed but different R making it smaller until the answer will not depend on R. 



We use the solution given by Eq. ( 2.21 ) , as a check of the accuracy of our numerical calculations. The 
procedure of the derivation of Eq. (2.33) from Eq. ( p. 15 ) is standard and it is described in Ref. [8] for 
example. 

Fig. |l| shows the numerical solutions of the BFKL equations (see Eq. ( |2.15D and Eq. ( |2.33| )) for 
d In (m bfkl (l)Y)) jdl (see Fig. 0-a) and In (m bfkl (I; Y) ) (see Fig. @-b) as functions of I = In J^^qcd 



for different values of the short distance regulator R at Y = 43. Comparing these curves with the prediction 
of the solution of Eq. (2.21) shown in Fig. || by the red curve, one can see that our numerical problems are 
concentrated in the region of very short distances where the value of x\2 approaches R. However, we see 
that in the region — 20 < / < 20 the numerical solution coincides with Eq. (2.21) with good accuracy. 
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Figure 2: dhvij\f BFKL (l;Y) )/dY (see Fig. ||-a and the zoomed picture in Fig. ^-b) for the BFKL equations 



(see Eq. ( 2.15 ) Eq. ( 2.33 ) with the BFKL kernels) as function of Y at I — 0. as is chosen to be equal to 0.2 and 
I = In {x\ 2 Aqcd)- Curves correspond to different values of R. The red curve shows the solution of Eq. ( 2.21 ). 



In Fig. g we plotted d In (j\f (/; Y) I jdY as a function of Y for different values of R. Comparing 

the curves for different R with the solution of Eq. ( |2.21| ) shown by red line in Fig. |2|one can conclude: first, 
that the value of the BFKL Pomeron intercept does not depend on the value of R and, second, that the 
numerical solution reproduces quite well the analytical one, given by Eq. ( 2.21 ) . However, Fig. g-b shows 
that our numerical value for the Pomeron intercept systematically above the analytical one approximately 
on 0.001 which we consider as a systematic error of our calculations. 



3. Modified BFKL Pomeron 



In this section we solve Eq. ( |2.15| ) and Eq. ( p. 16 ) in which the BFKL kernel are replaced by K (X13, x^2\x\2)- 
As has been mentioned above the numerical calculations in Refs. [14-17] show that such a modification of 
the BFKL kernel leads the exponential decrease of the scattering amplitude at large values of the impact 
parameter. Actually, we can see this directly from the equation (see Eq. (|2.14 )) . Indeed, one can see that 
the main contribution at large b stems from the region where \b — x%2\ < ^12- At such X32 the equation 
takes the form 

dN (X12 ra 26, b;Y) ._/",<, 1 



dY 



7 



Aa s / d X13 



-4Bb 2 -Bx? 



46 2 



2AT(xi 3 ,0;y) ex e 



-4Bb 2 



(3.1) 



We can replace exp (— B x\%) by Kq (/XX13) to reproduce correct exp (— /u6) at large b (see more detailed 
analysis of the form of the BFKL kernel in Ref. [14]). Nevertheless, we solve the equation with the kernel 



of Eq. ( |1.2[ ) since we interested in the behaviour of the (\b 2 (Y, I) |) versus energy for which the particular 
form of the kernel is not important. 

The initial condition for solving Eq. ( p. 15 ) and Eq. ( 2.16| ) we take in the following form 



N BFKL (x 12 ; Y = 0) = In (1/ {x\ 2 B)) ; Kf BFKL (x 12 ; Y = 0) = b 2 In (1/ [x\ 2 B)) ; (3.2) 

with b 2 Q = 1/B. 



It should be stressed that Eq. (2.18Q follows from the calculation of the amplitude for dipole- dipole 



scattering calculated in the Born Approximation of perturbative QCD. Introducing cutoff B we still consider 
the same initial conditions, that follows from the perturbative QCD calculation, since our main goal to 
study the influence of the evolution in Y with the modified kernel of Eq. (|1.2|) on the behaviour of the 
scattering amplitude at large b. 

It should be stressed that using the initial conditions of Eq. ( j3.2| ) one can see that introducing xf k = 
B xf k the modified kernel takes the form 

K B (xi3,x 32 |xi2) — ► K 1 (X13,X 32 |X12) (3.3) 

while the initial conditions of Eq. (|3.2| ) can be re-written as J\f BFKL (x 12 ; Y = 0) = In (l/ (xf 2 )j- 

Therefore, in xtk the equation and the initial conditions have the same form for any values of B. 
However, we prefer to solve equations with the kernel K B (xi3,X3 2 |xi 2 ) using the B independence of the 
solution as the way to check the accuracy of our calculations. For numerical solution we again introduce 
the short distance regulator in the same way as in Eq. ( p. 33 ) replacing K B (£i3,X3 2 |xi 2 ) by the following 
expression 

Jd 2 x 13 K B (x 13 ,x 32 \x 12 )Af BFKL (x 12 ;Y) = (3.4) 

I d 2 x 13 e ~ B ^f [2N BFKL {x 13 -Y) - 2 2 i\ ,„ 2 M BFKL {x 12 ;Y)\ 

J £32 + R z I xf 3 + x 23 + 2 R z ) 

3.1 Pomeron intercept 

3.1.1 Numerical solution for the Pomeron inrtercept 



In Fig. H it shown the solution to Eq. ( p. 15 ) with K^ (xi3,X3 2 |xi 2 ) for different values of B. One can see 



that N {x\ 2 ;Y) grows as exp(o'ol r ) and the solution does not depend on the value of B. In this figure 
we do not show the dependence on the value of the regulator R but we actually studied this dependence 
in the same way as for solution to the BFKL equation and saw that the solution does not depend on 
the value of R . To see the power-like dependence of the solution in a clearer way we plot in Fig. || 
din (J\f (xi 2 ; Y)) jdY . We see that the value of the intercept is smaller than the BFKL intercept but it is 
still increasing approaching this value. In Fig. [| we show that I dependence of the intercept is very similar 
to the BFKL Pomeron. 
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Log[N[Y,0]] 




Figure 3: The solutions to Eq. ( |2.15D 
with Kft (xi3,x 3 2|xi 2 ) ior Af (xi 2 ;Y) as 
function of Y. a$ is chosen to be equal 



to 0.2. I = In I 
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equation with the BFKL kernel (B = 



in the figure) / = In [x\ 2 ^ 
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N[Y, 0] 
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B=0GeV 2 


0.548 








0.546 







B=lGeV 2 


0.544 








0.542 







Fig. @-b 



Figure 4: d In J\f (xi2'> Y) /dY for the solutions to Eq. ( 2.15 ) with K^~ x {x\z, £321^12) as function of Y, as is 
chosen to be equal to 0.2. I = In (x\ 2 B) = 0. For the equation with the BFKL kernel (B = in the figure) 
I = In {x\ 2 A 2 QCD ) = 



3.1.2 Variational method 

The numerical solution suggests that the intercept of the Pomeron at B 7^ is the same as for the BFKL 
Pomeron. This fact as well as the influence of cutoff B on ujq can be understood from the BFKL kernel in 
7-representation (see Eq. (|2.2[) ). Indeed, this kernel appears in the calculation as the following integral ( 
seeRefs. [8,10,11]) 



X(7,0) 



1 00 

£7-1 

dt h 

1-t 



dt- 



t*i- 



t-1 







|t-l| V4t 2 + 1 



(3.5) 



where t = x\^/x\ 2 - 

The first term describes the increase of the dipole size due to evolution while the second one corresponds 
to normal DGLAP evolution in which the dipoles sizes decreases with the growth of Y. Introducing the 
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1 



#[38, 1] 



<9yN[38,l] 




— B=0GeV 2 

— B=2 GeV 2 

— B=4 GeV 2 

— B=8 GeV 2 



Figure 5: d InAf (x 12 ; Y) JdY 

for the solution Eq. (2.1 5|) with 



Kr { x i3: x 32\xi2) as function of I. 
as is chosen to be equal to 0.2. 
I = In {x\ 2 B) . For the equation with 
the BFKL kernel (B = in the figure) 
l = ln(xl 2 A 2 QCD )=0. F = 38. 



modified BFKL kernel we cut the sizes of the intermediate dipoles such that x 2 3 < 1 (t < l/x\ 2 . Assuming 
x\ 2 < 1 we see that we do not change integration over t < 1 but t should be smaller than t < l/x\ 2 . 
Therefore, we can estimate the modified kernel using 



X (7,0; £i 2 ) 



dt 



P 



-1 



l/«?a 



1-t 



dt- 



V 



-1 



t-1 



V^ia 1 

* 



1 



1 



|t-l| V4t 2 + 1 



1 

Introducing the variable t = x\ 2 jx\ 2 which is smaller that 1, we can re-write Eq. (3.6) in the form 



x(7,0;xi 2 ) 



, f<- 1 - 1 



dt 



t"T 



1 



+ 



rfi 



1-* A/x?a t\/4t 2 + l 



(3-6) 



(3.7) 



Taking the integral we obtain 



x(7,0;x 



12 



.BFKL 



X 



(7, 0) + arccsch {2/x\ 2 ) - 5 (xi 2 ; 1 - 7, 0) - In (l - x\ 2 ) 



(3.8) 



where 5 (£12; 1 — 7) 0) is incomplete .B function and ^ BFKL (-^ 0) is given by Eq. Q2.2J ). 

One can see that Eq. (3J3) depends on x\ 2 through £12 and, therefore, ujq (7,0,2:12) = asx{l-,^ x i2) 
cannot be the intercept of the Pomeron since the intercept cannot depend on the sizes of dipoles. However, 
we can use Eq. Q3.8| ) for developing the estimate in the variational approach for finding the ground state 
(the maximal intercept). Indeed, if we introduce^ 



N 



BFKL 



(x 12 ;Y) = |^e^*(x 12 , W ) 



the equation looks as 



CJ^(x l2 ,U)) = d 2 Xi 3 K B (xi 3 ,X 2 3\xi 2 )^ (Xl 3 , UJ) 



(3.9) 
(3.10) 



3 In this section as well as in sections 3.1.3 and 3.1.4 we use variable 3^ which is equal to asY(y — qjF). 
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where the r.h.s. of the equation has been discussed in Eq. ( |3.4| ), For the BFKL equation the eigenfunction 
v £ bfkl (#13) w) = (^13) an d the eigenvalue is given by Eq. ( |2.2[ ). This function cannot be a eigenfunction 
of Eq. ([T^) at -B 7^ as we have discussed. However, we can use this BFKL eigenfunction as a trial function 
in the variational method for searching the maximal value of the intercept. Actually, we use as the trial 
function 



^12) 



*tp(5l2,7) = ® BFKL (X13, 7) ®(1-^12J = (^13) ©i 1 

The variational principle has the following form for the problem of finding the maximal intercept 

/ d 2 x 12 ^* tr {xi2,"i) J d 2 xi 3 K B (»i3, £23^12) *tr (£13, 7) 

w m az ^ w (7) — r ,9- — , Tl * ,~ r~7. — jz r 

J d 2 X 12 *tr (Xi2,7) ®tr (^12,7) 

^ f rf 2 Xi 2 ^* r (gig, 7) X (7, 0; S12) g«r (^12, 7) 
Jd 2 Xi2^ t * r (xi2,7) *tr(»12,7) 

From convergency of J d 2 x\2 *| r (£12,7) ^tr (£12,7) we conclude that 7 > 0.5. 




Figure 6: w (7) (see Eq. (3.12)) and 
omBFiCL (7)(see Eq. (|2~2])) versus 7( 
black and red lines, respectively). In 
this figure we redefine w dividing it by 

as- 



(3.11) 



(3.12) 



From Fig. y we see that w(7 = 1/2) coincide with the BFKL value. In other words, we prove that the 
resulting maximal oj can be either equal to the BFKL one or larger than the BFKL value. 

3.1.3 Semi-classical solution 

Examining the property of the solution to the modified BFKL equation we wish to find a semi-classical 
solution to this equation searching it in the form: 



N(y-l) = e s ^ = e "W)Y + (7W)-W whe re u>(y,l) = 9S ^ ° ; j(y,l) 



oy 



1 



dS\(y;l) 
dl 



(3.13) 



with smooth functions oj (y, I) and 7 (y, I). 

Inserting Eq. ( 3.13| ) to the equation we obtain: 

00 (y, I) - x(7,0,x 2 2 ) 







(3.14) 



Deriving Eq. ( |3.14 ) we use that 7 is a smooth function and performing the integral of Eq. j p.Y|) we can 
consider 7 as being a constant. Since the equation is the first order differential equation in respect to Y 
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the condition for applying the semi-classical approach looks as follows 

^m^L « 1; z ^m « ! (3 . 15 ) 

(l- 7 (^,/)) 2 dl 

where/ = \n{x\ 2 / Kq C d)- 

It is known [36] that for the equation in the form 

F(y,l,S,-y,u>) = (3.16) 

we can introduce the set of characteristic lines : l(t),y(t),S(t), uj(t) and "f(t) which are the functions of 
the variable t ( artificial time), that satisfy the following equations: 

{ '' dt 7 dj 



dS d X (7, 0,0 

-dl = 'T E r + uF - = -(T- 1 )— ^~ 



(4)g = -( fl+7 F,) = MrWj!MW) (3 . 17) 



The fifth equation is Eq. ( pl4| ). 

First, we see two thing which simplify a bit the equation: we can consider t = y and dividing Eq. ( |3.17| )- 
4 by Eq. (|3.17|)-1 we see that 



d 7 _ m 



(3.18) 



^ dx(-y,0,l(t)) 

d-y 

We solve this equation putting j(l = oo) which coincide with the solution to the BFKL , as the initial 
condition. Due to convergency of the integral for the norm of N we know that 7 (00) > 0.5. 

The second step after finding 7 (I) is to solve Eq. ( |3.17| )-1 to find the form of trajectories. It tuns out 
that we do not need to find I (y) as a function of y for finding the values of oo. 

The last the third step is to find from Eq. ( 3.14 ) the value of u. In Fig. we plot function 7 (/) for 



different 700 = 7 (/ = 00), while Fig. g shows the dependence of uj (7) = uj (y, I) — x (7> 0; versus I. 

From Fig. || one can see that the values of 10 (7 (I)) does not depend on I and these values turn out to 
be smaller that for the BFKL equation. 

In Fig. H we plot the ratio d'y (I) /dl J (1 — 7 (/)) and the product Idj (I) /dl as a function of I. One 
can see that both of these observable turn out to be small and, therefore, we can trust the semi-classical 
approach. 

Concluding this subsection we see that the semi-classical approach leads tow = ojbfkl- 
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700=0.6 
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Figure 7: The trajectories of Eq. (3.14) at different values of 7(00) 
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Figure 8: u (I) versus I at different values of 7 (00). Notice that w B pkl (0-5) (blue line) and to (700 = 0.5) (red dotted 
line) are the same. 



3.1.4 Diffusion approximation 

In direct analogy with the BFKL equation we develop in this section the diffusion approximation to the 
modified BFKL equation. The main idea of this approximation is to introduce a new function 



N(y,l) = e~2 l N(y,l) 



e+ioo J. . 
e-ioo z,/u ' 



e+ioo 



dui 



le+oo 



du 



(■—100 27T2 Jie— 00 ^ 7ri 



-n [u, v) e 



u)"y + iul 



(3.19) 



The diffusion approximation means that we can reduce the modified BFKL equation to the differential 
equation using the following expansion for n(uj,l): 



1 ,/\ / ,\ dn(oj,l'), ,,, A ld 2 n(oj,l') . ,-, ,\2 



dl' 



2 dl' 2 



(3.20) 
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Fig. |-6 



Fig. g-b 



Figure 9: The ratio d"f(l)/dl/ (1 —7(0) (see Fig. ^-a) and the product ldj(l) /dl (see Fig. ^|-b) versus I for 
7 (oo) — 0.5. 



The equation takes the form after plugging in Eq. ( |3.2C| ) into Eq. ( p. 10 ) 



/ a /,m / ,\ ,-, ,,<, dn(co,l) 1 ,„,,. d 2 n(uj,l) 

(u - A(0)n(w,0 - dl(i) ^ ; - -d2(i) d K p ' 







(3.21) 



Functions A (Z),dl (/) and d2 (I) can be expressed through the kernel % ( | + ^ v i 0\ e2 ) that has been 
introduced in Eq. (|3.8|) , viz, 



Ai/i - v (l,0,e^): r/J.(/): /- 



cV 



-' l d 2 x [\ + iv,0,e 7 < 

V=Q - d2(l) = ^— '- \ v=0 . (3.22) 



c>z/ 



Fig. 1C shows these functions. 

All) 



did) 



(12(1) 




Figure 10: Functions A (Z), dl (Z) and <i2 (/) versus /. The red line shows the BFKL value for w BFK l (I) = 41n2. 



Introducing n (to, I) = exp (eft (u, I)) for function (ft Eq. ( |3.21| ) reduces to 

(w-A(O) = dl(07(w,0 + |rf2(0 (^|^ + 7 2 (^,0 



(3.23) 
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where 7 = ^' ' . 

For I -> -00 A (I) -> w BFKL , dl (Z) ->■ and d2 (I) ->• L> (see Eq. (|2T2C|)). 
Therefore, for large and negative I the solution for 



7(w,0 



vT 



W — CJt: 



/(2D ). 



These values give us the initial condition for Eq. ( 3.23| ). 



Solving Eq. (3.23) numerically we see that for all u < w B fkl 7 < at Z 
amplitude n which normalization is convergent integral, namely, 



(3.24) 



+00. Such 7 lead to the 



/ 



dl\n (u), I) \ 2 dl < 00 



(3.25) 



while for uj > w B pkl 7 at large I > is positive and the integral in Eq. ( |3.25| ) is divergent (see Fig. |j~l| ). 
In other words, our spectrum of u) is continuous with u) < o; BFKL . 
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Figure 11: Solution to Eq. (3.23) with the initial condition of Eq. (3.24) versus I: Fig. [Tlj-a for Re7(w,Z) and 
Fig. ||-bforlm7(w,0. 



3.2 Pomeron slope 



Approaching the problem of calculation of (|6 2 |), the first observation that we can make is that Eq. (2.16) is 
valid also for the modified BFKL kernel. Indeed, the kernel itself does not depend on the impact parameter 
and functions E n '^ (pi a , P2a) is the complete set of functions. As we noticed in derivation of Eq. fl2.16| ) 
the term in < . . . > in Eq. ( |2.17| ) vanishes due to invariance of function _E n=0 >7 (pi a ,p2 a ) with respect to 

transformation b — > — b. Since uj (7, ri) < for n > 1 the only component of the arbitrary function of 
the initial condition that survives at large Y is its projection on E n=0,,y (pi a ,p2a) ■ F° r this projection the 
term of Eq. ( 2.17 ) vanishes leading to Eq. ( p. 16 ). 
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Figure 12: {\b 2 (Y, I) |) versus Y at different values of B (see Eq. (jT|)). a s = 0.2. For both B 
we took &o in Eq. (2.24) and Eq. (3.2) the same and equal to 0.5 GeV~ 2 . 



and B = IGeV 2 



Solving Eq. ( 2.16| ) with the kernel of Eq. (1.2) we calculate (see Eq. ( 2.1 3| ) ) (|& 2 |) as function of Y (see 



Fig. |12j). In the Reggeon approach this value gives the information on the slope of the Pomeron trajectory 
(a'jp) since (|6 2 |) > 4 a'jp Y . One can see from Fig. 12 that {\b 2 \) does not depend on Y at large values 



of Y. Actually, the solution to the modified BFKL equations shows a weak Y dependence (see Fig. 13 which 
is zooomed Fig. 12). However, even if we assume that (|6 2 |)|^ =1 = Aa'pY the value of a'jp < 0.510~ 5 &o 
is extremely small. Hence we can conclude that the modified BFKL Pomeron has a'jp ~ 0. 
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Figure 13: (|6 2 (Y, I) |) = Aa'p Y versus Y at different values of B (see Eq. fll.Ifl). as = 0.2 



This result was expected and its explanation based on the general features of QCD. The general origin 
of the increase of (\b 2 \) with Y ( The relation: (|6 2 |) = Aa'jpY) , was understood in 70's by V.N.Gribov 
(Gribov's diffusion [40]). Each emission leads to change in impact parameter by A6 2 = 1/p 2 - where px 
is the typical transverse momentum. After n emission (6 2 ) = A6 2 n which corresponds to the random 
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Y=50, B=lGeV 



Y=200,B=lGeV 2 



-60 -50 -40 -30 




Figure 14: (\b 2 (Y,l)\) versus I = 
In (x\ 2 B) at different values of Y, Solid line 
describes Y — 50 while the dotted one cor- 
responds to Y = 200. a s = 0.2. 



walk in the transverse plane. Since the number of emission is proportional to Y we obtain (|6 2 |) oc Y. 
In the parton model the typical pr is independent of Y and, therefore, we see that diffusion in b leads to 
«jp 7^ 0. However, in QCD average pt depends on Y. Such dependence stems from the diffusion in Inpx 
which means that in each emission of gluons In p? changes by a constant. Being a general features of all 
theories with dimensionless coupling such diffusion in liipx comes out from the BFKL equation leading to 
(In (pt/Pq x) = ^ Ad Y ( see Eq. ( 2.20Q ). From this formula one can see that we see two different branches 
in the BFKL equation: one leads to a rapid increase of the typical transverse momentum while another to 
a steep decrease. This decrease does not influence the calculation of the average pt, resulting in a'jp = 
for the BFKL equation. Modeling confinement by introducing cutoff in the sizes of produced dipoles we 
prohibit the decrease of the typical transverse momenta of the emitted gluon. As a result, the only diffusion 
in the large transverse momenta occurs leading to negligible a'jp. 



Fig. 14 shows the dependence of (\b 2 (Y,l) |) on / = In (x\ 2 B\ at B = lGeV 2 . We can see that the 
dependence of (\b 2 (Y,l) |) on I is rather weak except the region of I close to 0. The maximum at / = 
stems from Gribov's diffusion during the first several emissions until the average pt grows to a considerable 
value. 



3.3 Saturation momentum 

We have demonstrated that the modified Pomeron has a correct behaviour at large b but violates the 
s-channel unitarity (Froissart theorem [9]) both for the partial amplitudes and for the total cross section, 
since they are proportional to s u ' BFKL . Therefore, we need to develop the CGC/saturation approach [5-8] 
and reference therein), based on the modified BFKL Pomeron to obtain the amplitude that will satisfy 
the unitarity constraints. We are going to develop such an approach but in this paper we wish to use the 
well known feature of the CGC/saturation approach: the energy behaviour of the new dimensional scale 
(saturation moment) can be found from the linear equation (see Refs. [5,37,38]). This scale is the solution 
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■,[Y] 




Figure 15: l 8 (Y) = In (q 2 s (Y) / Q 2 a (Y ) 

, where Q s (Y) is the solution to Eq. ( 3.26J ), 



versus Y for different values of Ao- as = 0.2. 
The red solid and dotted lines correspond 
to solution of the BFKL equation, the solid 



dark brown line describes Eq. ( 3.27 ) while 
the blue line is the solution to Eq. ( 3.26 ) for 
the modified BFKL equation. 



of the equation 



Af- 



BFKL 



\Y 



■A/n < 1 where 



An = Const 



.Qs(Y)> 
For the BFKL equation the solution to Eq. (3.26) is known. It takes the form 



l s (Y) = In (Q 2 (Y)/Q 2 (Y q ) 

Uj(j C r) „ r ... 3 



(3.26) 



(3.27) 



l-7c 



(X ~ Y ) 



2(1 - 7c 



HY/Yo 



2vr 



(1 - 7 cr) 2 V W "(7cr) 



Y 



Ta 



where Y = ln(l/x) is our energy variable, oj"(j) = d 2 Lu(-y)/(d"f) 2 , the value of j cr can be found from the 
equation [5,37]: 



w(7 cr ) 



d^cr 



, with cj (7) = asx (7) = as (2ip (1) - ip (7) - ip (1 - 7)) 



(3.28) 



where ip (7) = dlnT (7) /drj and T (7) is the Euler gamma function. In Eq. ( p.27| ) the first term was found 
in Ref . [5] , the second in Ref . [37] and the third term was calculated in Ref . [38] . The solutions to Eq. ( |3.26D 
for different values of Ao are plotted in Fig. [D| Two features are clear from Fig. 15: Eq. (|3.27f) is in a 
good agreement with the numerical solutions; and the energy dependences of the saturation scale are quite 
different for the BFKL equation and for the modified BFKL equation which includes confinement. 

One can see that the energy dependence of the saturation scale for the BFKL Pomeron with confinement 
(B = 1 in Fig. ^) turns out to be much milder than for the BFKL evolution (B = in Fig. 15). 



4. Conclusions 



The main goal of this paper is to find how our assumption that the size of produced dipoles cannot be 
large, will affect the main properties of the BFKL Pomeron. To achieve this goal we solved the BFKL 
equation with the modified kernel of Eq. ( |l.2|) . We found out that the modified BFKL Pomeron has the 
same intercept A as the BFKL Pomeron ( Agpj^ = uj bfkl = 4 ln2a,g) and a'jp = 0. Therefore, the 



20 



BFKL Pomeron with the modified kernel reproduces the main features of the soft Pomeron that has been 
found both from N=4 SYM theory [32-35] and from the high energy Reggeon phenomenology [29,30]: the 
large value of the Pomeron intercept (loq ~ 0.2 — 0.3) and a'jp = 0. These both conclusions are in agreement 
with the numerical solution of the modified BFKL and BK equations [14,15]. 

We consider as one of the results of this paper that we developed several methods to solve the modified 
BFKL equation analytically (semi-classical and diffusion approximations). The fact that these methods 
work we checked with the numerical calculation. 

Actually, we were surprised that the model for confinement changed so little in the BFKL Pomeron 
and on qualitative level, the Pomeron that emerges from the modified BFKL equation, looks quite the 
same at the BFKL Pomeron, both in parameters and in character of the energy behaviour. It seems that 
the only difference between the BFKL Pomeron and the modified BFKL Pomeron is that the second has 
a correct large impact parameter behaviour. The most sensitive to the confinement global feature of the 
BFKL Pomeron turns out to be the energy behaviour of the saturation momentum. The confinement 
diminishes the energy dependence of Q S (Y) bringing it closer to the experimental one. 

Eq. ( 2. IE ) is new and using this equation we are able to calculate directly {\b 2 (Y, l)\) which gives the 



information on the effective slope of the resulting Pomeron. 

We believe that this paper will be useful in the search of the theoretical motivated way to include the 
non-perturbative corrections at large values of the impact parameters as well as in understanding of the 
main ingredients of high energy phenomenology for soft processes. In the future publication we are going 
to study how the way of introducing confinement into the BFKL equation could change the features of the 
Pomeron and to develop the CGC/saturation approach, based on the modified BFKL Pomeron. 
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